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Abstract 

The decay of colour neutral 2-gluon systems in qq and 2q2q has been simulated 
with the Monte-Carlo method, taking into account an effective 1-gluon exchange 
interaction between the emitted quarks, which was folded with a 2-gluon density 
determined self-consistently. Finite 2-gluon densities were obtained with a mean 
square radius < > of about 0.5 fm^. These give a significant contribution to the 
2-gluon field correlators and the gluon propagator from lattice QCD. 
By solving a relativistic Schrodinger equation the binding potential of the 2-gluon 
system is computed. By the relation 2g — > {qq)^ this corresponds to the confine- 
ment potential between quarks, which is in excellent agreement with the confinement 
potential from lattice QCD. The lowest bound state in this potential allows an in- 
terpretation of the 2-gluon mass as binding energy of the two gluons. The lowest 
glueball with Eg = 0.68 it 0.10 GeV is consistent with the scalar cj(550), with a 
mass significantly lower than obtained in lattice QCD simulations. However, the 
first and second radial excitations agree with the lattice results. 



PACS/ keywords: 12.38.Aw, 12.38.Lg, 12.39.Mk, 02.70.Uu/ gluon-gluon cou- 
pling in QCD, eff'ective 1-gluon exchange interaction, lattice gluon field correlators 
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Quantum chromodynamics (QCD) is not well understood in the infrared region, where 
the most important properties of the strong interaction are generated, i.e. the hadronic 
masses and the confinement of quarks and gluons to hadrons. As perturbative methods 
are generally not applicable, simulations of the QCD equations on the lattice [ij together 
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with studies of Dyson-Schwinger equations [2, 3j offer the best known methods to study 
the structure of QCD in the infrared region. However, because of their complexity and 
severe problems in the lattice simulations (due to discretisation, with small quark masses 
and the need of a very high computing power) it is very difficult to reveal the underlying 
structure of QCD. Therefore, alternative approaches are welcome, which can shed further 
light on the underlying physics. 

In the present paper we discuss a phenomenological description of two-gluon systems 
following the old idea (see e.g. in ref. [I]), that the non-Abelian structure of QCD may 
allow a coupling of two gluons to bound colour singlets with spin J=0 (and 2). For the 
description of such a scalar bound state $ we may write the radial wave functions in the 
form 



where ipj{rj) are the radial wave functions of the two gluons. To investigate the properties 
of such a 2-gluon system we studied the decay 2g — » qq+2q2q with an attractive interaction 
between the emitted quarks. 

There is evidence for the existence of scalar and tensor gluonic fields from the formation of 
glueballs in pure Yang-Mills theory [5] , from the mechanism of Pomeron-exchange P, [7] 
in high energy hadron-hadron scattering (which is understood by the exchange of two 
non-perturbative gluons which couple mainly to spin zero), scalar excitation of baryon 
resonances j8j and hadron compressibilities Finally, the effective confining quark-quark 
potential fTO] contains a sizeable scalar contribution. 




Assuming an effective 1-gluon exchange interaction Vig{R) = —ag/R between the emitted 
quarks (or antiquarks) with relative distance R = |rj— rj|, the decay from a 2-gluon system 
2g —>■ (qq)^ requires a modification of the free q-q interaction by the density of the 2-gluon 
system, which may be expressed by a folding integral 



where p$(r ) is the 2-gluon density p$(r ) = |'?/'$(r )p. 

It is interesting to note, that for a spherical density the Fourier transform of eq. ([2]) to 
momentum (Q) space yields 
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where p<s>{Q) = An J r'^dr jo{Qr) p<i>{r). Comparing this with the standard 1-gluon ex- 
change force yields a Q-dependent strong couphng 

as{Q) = as p^{Q) , (4) 

which is quahtatively consistent with the known fact of a "running" of as{Q) and the 
condition as{Q)-^ for Q^oo (asymptotic freedom) 3- Note, that the 2-gluon densities 
discussed below yield as{Q) in reasonable agreement with the lattice data of ref. [12j for 
small momentum transfers. 

Further, the interaction has to be modified by the finite size of the decaying 2-gluon sys- 
tem. This is approximated in our calculation by replacing the 1-gluon exchange potential 
Vig{R) by 

Vl.iR) = V,,{R) e^-'^^^) (5) 

with the requirement < r^, >^< r^^ >. This effect is shown in the upper part of fig. 1, 
which gives a 2-gluon density by the upper dot-dashed line and pure 1-gluon and modified 
1-gluon exchange potential by the lower dashed and dot-dashed lines, respectively. The 
resulting interaction is given by 

= / drp^{r ) Vi^{R - f ) , (6) 

which has the characteristcs of a quite rapid fall-off to zero for large radii (solid line in 
the upper part of fig. 1). 

Finally we have to take into account the fact, that the decay into 2g2g favours relative an- 
gular momentum L=0 between the emitted quarks, whereas for the decay in qq the outgo- 
ing quarks are in a relative p-state (L=l). For the p-wave density p%{r) = p%{r) Yi^rn{6, 4>) 
we take p$(r) = (1 + /? ■ d/dr)p^{r) with the constraint < r > = j dr rp^ecij) = (to 
suppress spurious motion), by which (3 is fixed. 

By relativistic Fourier transformation [13] the effective interaction (I6l) can be transformed 
to momentum space 

= 4^ / R'dR Jo{QR) V^^iR) , (7) 

with Q' = Q \Jl + [Q"^ /4mW and m$ being the mass of the 2-gluon field. For p$(r) and 
p%{r) the potential in momentum space V^q(<5) is given in the lower part of fig. 1 by the 



^However, is is clear that the 2-gluon densities deduced here fall off much faster than the slow loga- 
rithmic decrease found for as{Q) (see [llj). 
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solid and dot-dashed line, respectively, which show significant differences. Due to the 
finite size effects both interactions are finite at vanishing momentum transfer. 

Using this effective q-q interaction Monte-Carlo simulations of gluon-gluon scattering have 
been performed in fully relativistic kinematics using the CERN routine GENBOD [11], in 
which the 2 gluons in the final state can decay in qq and 2g2g (we used massless quarks). 
Random Q-transfers between and 6 GeV/c have been used, from these the momenta Pi 
of the outgoing quarks or antiquarks were also determined randomly, but restricted by 
the available phase space. The potential Vgg(Ap) ([7]) has been used as a weight function 
between the outgoing quarks (with the relative momenta Ap = pi —pj) in addition to the 
normal phase space weight. Resulting gluon momentum distributions dqq{Q) and d2q2q{Q) 
for decay into qq and 2g2g were generated. Their sum D^{Q') = dqq{Q') + d2q2q{Q') can 
be related to the radial density p$(r) of the 2-gluon system 



with Q' as in eq. (171). The mass m$ to determine Q' in the relativistic Fourier transforma- 
tions ([7j) and ([8]) has been used as a fit parameter, for which a value of 0.68 GeV yields 
the best results. 

The condition, that p$(r) in the interaction ([6]) and in eq. ([8]) should be the same, al- 
lows us to determine this density. As the 2-gluon density was not known initially, our 
simulation procedure had to be iterated: starting with a certain density p$(r) the folding 
potential ([6]) was calculated. The Fourier transform ([7]) was then inserted in the Monte- 
Carlo simulations, from which gluon momentum distributions dqq{Q) and d2q2q{Q) were 
generated as shown in the upper part of fig. 2. By inversion of eq. (IHD a Fourier retrans- 
formation of Dq,{Q) to r-space had to be made and the resulting density p$(r) compared 
to the one used in eq. ([6]). Self-consistency of both densities was required. 

Resulting momentum distributions dqq{Q) and d2q2q{Q) for a self-consistent solution (dis- 
cussed below) are given in the upper part of fig. 2. These are very different for the two 
decays: the decay in qq falls off rapidly with increasing Q, whereas a much slower fall-off 
for 2q2q decay is observed. This difference can be readily explained by the fact, that 4 
quarks can absorb much higher momenta than 2 quarks. We see that the sum D^{Q) is 
in reasonable agreement with p<s>{Q) from the Fourier transformation of p$(r) inserted in 




(8) 
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eq. ([6]) (dot-dashed line), which is quite well approximated by a radial dependence 

^<s>ir) = exp[-{r/aY] (9) 

with values of a and k in table 1 and ipo taken from the normalisation 47r / r'^dr ^/'$(r)^=l. 

The self-consistency condition for our simulations is far from trivial. The resulting density 
is given in the lower part of fig. 2, which indicates clearly that a self-stabilized 2-gluon 
field is generated. In table 1 we give two solutions of about equal quality, which may 
reflect the minor ambiguities of our method. The dashed area indicates the estimated 
uncertainties of the deduced density. 

The extracted 2-gluon density should give a significant contribution to the gluon 2-point 
functions extracted from lattice QCD simulations in form of gluon field correlators p!5l[T6] 
and the QCD gluon propagator (see e.g. [IZj). This is shown in fig. 3. In the upper part we 
make a comparison of our results with the 2-gluon field correlator C_[_{r) of Di Giacomo et 
al. [16]. C^{r) has been weighted exponentially with distance (yielding < >~0.01 fm^); 
therefore we have to compare the logarithmic values of C±{r) with our 2-gluon density 
(lower solid line). This leads to a reasonable agreement at large radii. For a comparison 
with the gluon propagator we have taken the lattice data (obtained in Landau gauge) of 
Bowman et al. [T8]. These are compared with the Fourier transformed 2-gluon density 
p$(Q) shown by the lower solid line in the lower part of fig. 3. 

In spite of a reasonable agreement with our densities at large radii or small Q-values, 
respectively (note that in fig. 3 the gluon propagator is multiplied by Q^), both lattice 
data show deviations, which may be described by a 2-gluon vector component (where 
the two gluons couple to angular momentum L=l). Such a (current) field is constrained 
by the condition < r >= to suppress spurious motion. We assume a form pvec{f) = 
Po(l — 0.25 r/a') exp{ — {r /a')} with a'=0.11 fm, which fuUfills the requirement < r >= 0, 
yielding a contribution given by the dot-dashed lines in fig. 3. Using for this a strength 
of about 18 % of the scalar component yields a consistent description of the (gauge 
independent) 2-gluon correlator and the gluon propagator in Landau gauge. The latter is 
not very different in Laplacian gauge [19], but a fit to the gluon propagator in Coulomb 
gauge [20] yields a relative vector component about a factor 2 smaller, but in this gauge 
the gluon propagator is not covariant. 
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Table 1: Deduced density parameters and energies of the lowest radial eigenmodes with 
the units fm, fm^ and GeV for a, <r^ >, 6m and Ei, respectively. 



2-gluon field 


K 


a 


<r2 > 


6m 




Eo 


El 




E2 


(solution 1) 
(solution 2) 


1.50 
1.51 


0.66 
0.69 


0.48 
0.52 


0.04 
0.12 


0.6< 
0.6< 


BiO.lO 
BiO.lO 


1.70±0.15 
1.74±0.15 


2.5< 
2.6< 


B±0.20 
B±0.20 



In various investigations of the gluon propagator and Pomeron-exchange it has been con- 
cluded, that 2-gluon fields should be massive. Cornwall [1] has discussed the gluon 
propagator in relation to a (relativistic) generation of mass. In studies of Pomeron- 
exchange [HI [7] an effective gluon mass between 0.3 and 0.8 GeV has been extracted. In 
our relativistic Fourier transformations we had to assume m$ ~ 0.68 GeV. This is consis- 
tent with the gluon pole mass of 0.64 ± 0.14 GeV deduced in ref. [21]. We shall see, that 
the extracted mass can be understood as binding energy of the 2-gluon system including 
relativistic mass corrections. 

A finite 2-gluon density as shown in fig. 2 may be interpreted as a quasi bound state 
of the two gluons (glueball). Therefore, from the 2-gluon density the binding potential 
of the 2-gluon system can be obtained by solving a three-dimensional reduction of the 
Bethe-Salpeter equation in form of a relativistic Schrodinger equation 

- [7^ + - Mr))Mr) = E^Mr) , (10) 

where '?/'$(r) is the 2-gluon wavefunction ([1]) and /i$ a relativistic mass parameter, which is 
related to the mass m$ of the 2-gluon system by /i<i> = ;| m<j, + 6m. In the non-relativistic 
limit 6m = and /i$ corresponds to the reduced mass /i = ^^^^ which gives fi = \ m^ 
for mi = m2 = m$/2. Using a 2-gluon wave function of the form of eq. yields the 
explicit form for the binding potential 

v<,{r) = ^ [^/-r-' [<'-r 1)]] + Eo . (11) 

To calculate the shape of the potential the used value of Eq (which gives only an absolute 
shift) is not important. With k, a, 6m in table 1, slightly different solutions of the binding 
potential were obtained, which are given by the dot-dashed and dashed lines in fig. 4. 
Because of the relation 2g — ^ {qq}^ this potential can also be considered as confinement 



potential between the emitted quarks. This is in surprising agreement with the 1/r + 
hnear form expected from potential models p^j and consistent with the lattice data of 
Bali et al. [22]. It is important to note, that our potential ( fTTI) reproduces the 1/r + linear 
form without any assumption on its distance behavior; this is entirely a consequence of 
the deduced radial form of the 2-gluon wave function, see eq. ([9]). 

Absolute values of binding energies Ei in the potential ffTOj) can be obtained by fitting 
the binding potential by a form Vfit{r) = —^ + br, and adding to the eigenvalues a 
constant 6E to give a g.s. energy Eo = eo + SE consistent with m$. This yields bound 
state energies for the ground state and the two lowest radial excitations given in table 1. 

By computing also eigenstates in the q-q potential (jS]) we find only one bound state with 
an energy in the order of -10 MeV. This very low binding energy indicates clearly, that 
the glueball states must have a large width, since T 1/Eo should be valid for the decay 
width. This is consistent with the general expectation of the width of glueball states (F > 
500 MeV, see e.g. ref. [23]). From these results we may conclude, that the glueball ground 
state with i?o=0.68±0.10 GeV and a large width may be identified with the scalar cr(550), 
which is most likely not a qq structure (see e.g. the review by Bugg [24]). 

Glueball masses have been deduced also from lattice simulations 0,125], in which a glueball 
mass below 1 GeV has not been found. However, in these simulations 0^^ glueball masses 
have been extracted at about 1.7 and 2.6 GeV, which correspond very nicely to the first 
and second radial excitation in table 1. Our evidence for a low lying glueball is supported 
by QCD sum rule estimates [26], which also require the existence of a low lying gluonium 
state below 1 GeV. 

Conclusion: We have presented a phenomenological method which yields a good un- 
derstanding of important properties of the non-perturbative structure of QCD. The basic 
assumption, that stable 2-gluon colour singlets exist, which decay into qq pairs, has al- 
lowed a self-consistent determination of the densities of these 2-gluon systems. These 
give important contributions to the gluon 2-point functions, field correlators and gluon 
propagator, and yield a two-gluon binding potential which can be identified with the 
confinement potential determined from potential models and lattice QCD. 

A detailed comparison of our approach with lattice QCD will be important. In the latter 
two-gluon fields are facilitated by loop integrations. By relating gluon loops to our 2-gluon 
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fields (recall that Wilson loops are gauge invariant as colour singlet 2-gluon systems) one 
can imagine, that Wilson's confinement picture (by evaluating closed gluon loops around 
separated quarks) has a correspondence to our interpretation of confinement as binding 
of 2-gluons. 

New information is obtained on the mass of glueballs. Consistent with pole analyses of 
the gluon propagator and QCD sum rules we obtain a low mass of the 2-gluon system 
(glueball ground state) in the order of 0.7 GeV. This is in the mass region of the scalar 
£7(550). Higher radial excitations are obtained at about 1.7 and 2.7 GeV, consistent with 
0+"'" glueball masses from lattice QCD. 

We thank P. Decowski, S. Krewald and N. Nikolaev for fruitful discussions and M. Dillig 
for valuable comments and a careful reading of the manuscript. 
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Figure 1: Upper part: Two-gluon density p$(r) (upper dot-dashed line) and potentials: 
1-gluon exchange potential (dashed line), modified by size effects (dot-dashed line) and 
folding potential given by solid line. Lower part: Fourier transform of the folding 
potential using p$(r) (solid line) and p$(r) (dot-dashed line), respectively. 
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Decay of 2-gluon system 
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Figure 2: Upper part: Resulting 2-gluon momentum distributions (multiplied by Q^) for 
decay in qq and 2^2^ and sum. Lower part: deduced 2-gluon density with estimated error 
band. The dot-dashed hues represent solution 1 in table 1. 



11 




1 2 3 4 5 6 



Q (GeV/c) 



Figure 3: Gluon field correlator log{C±) from ref. [I6] (upper part) and gluon propagator 
from ref. [18] (lower part) from lattice QCD simulations in comparison with our results. 
The lower solid lines correspond to the density p$(r) and its Fourier transform p<i>{Q), 
respectively, and the dot-dashed lines to an additional vector component with the same 
relative normalisation in the upper and lower part. The sum of both contributions yields 
a good description of both lattice data. 
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Figure 4: Two-gluon binding potential (|TT1) given by the dot-dashed and dashed hnes (with 
the respective parameters from table 1) in comparison with the confinement potential from 
lattice calculations [22j. 
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